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ABSTRACT 

We analyse the dark matter (DM) distribution in a « 10 12 Mq halo extracted from 
a simulation consistent with the concordance cosmology, where the physics regulating 
the transformation of gas into stars was allowed to change producing galaxies with 
different morphologies. Although the DM profiles get more concentrated as baryons 
are collected at the centre of the haloes compared to a pure dynamical run, the total 
baryonic mass alone is not enough to fully predict the reaction of the DM profile. We 
also note that baryons affect the DM distribution even outside the central regions. 
Those systems where the transformation of gas into stars is regulated by Supernova 
(SN) feedback, so that significant disc structures are able to form, are found to have 
more concentrated dark matter profiles than a galaxy which has efficiently transformed 
most of its baryons into stars at early times. The accretion of satellites is found to 
be associated with an expansion of the dark matter profiles, triggered by angular mo- 
mentum transfer from the incoming satellites. As the impact of SN feedback increases, 
the satellites get less massive and are even strongly disrupted before getting close to 
the main structure causing less angular momentum transfer. Our findings suggest that 
the response of the DM halo is driven by the history of assembly of baryons into a 
galaxy along their merger tree. 
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1 INTRODUCTION 

Dark matter haloes in Cold Dark Matter (CDM) scenarios 
have been extensively studied during the last decades (e.g. 
Merrit et al. 2006; Prada et al. 2006; Gao et al. 2008; Navarro 
et al. 2008) , since understanding their properties is relevant 
not only for the study of galaxy formation but also for con- 
straining the nature of DM itself (e.g. Springel et al. 2008). 
Several controversies have been issues of thoughtful discus- 
sions such as the universality of the DM profiles (Navarro, 
Frenk & White 1996, NFW) or the systematically increas- 
ing central DM density with increasing numerical resolution 
(Moore et al. 1999). However, new improvements in the nu- 
merical codes and computational power have challenged pre- 
vious results. Recently, Navarro et al. (2008) found evidence 
for the non-universality of the DM profiles and for non-cuspy 
inner profiles by using high resolution N-body simulations. 

The assembly of baryons in galaxies is expected to drive 
changes in the DM distribution as it reacts to gas infall and 
outflows (e.g. Read & Gilmore 2005) and to the angular mo- 
mentum redistribution caused by infalling satellites (Barnes 
& Hernquist 1996; Debattista et al. 2008, D08). Blumenthal 



et al. (1986, B86; see also Young 1980 for a different imple- 
mentation) proposed the adiabatic contraction (AC) model 
to globally predict the effects of central baryonic accumula- 
tion on the DM distribution. Although it is clear that the 
condensation of baryons increases the DM concentration, 
it remains to be determined how this would proceed and 
how the transformation of baryons into stars and supernovae 
(SN) feedback might affect the DM distribution. Gnedin & 
Zhao (2002) studied the effects of maximum feedback on the 
central density of DM haloes of gas-rich dwarfs, finding them 
to be too weak to explain the observed rotation curves. In- 
terestingly, Read & Gilmore (2005) showed that two impul- 
sive mass-loss phases are needed to significantly change the 
DM profiles of dwarf galaxies, suggesting the need for bursty 
star formation histories. Also, several works have studied the 
change in shape of DM haloes toward rounder configurations 
induced by the presence of a central baryonic concentration 
(Kazantsidis et al. 2004). D08 claimed that a larger, irre- 
versible change in the DM distribution can be produced if 
there is transport of angular momentum from baryons to 
the DM. In a fully cosmological context, simulations have al- 
ready shown how the DM haloes concentrate with increasing 
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baryonic condensation, providing hints of a possible depen- 
dence on the assembly history of baryons (e.g. Tissera & 
Dominguez-Tenreiro 1998; Gnedin et al. 2004; Onorbe et al. 
2007). However, the high complexity of the problem requires 
high numerical resolution to describe in more detail the joint 
evolution of the DM and the baryonic distributions to get 
to more robust results. Also, none of these previous stud- 
ies included treatments for SN feedback and a multiphase 
interstellar medium (ISM) self-consistently. This is of par- 
ticular relevance since feedback not only regulates the star 
formation in the main galaxy, but also in the merging sub- 
structures and, as a consequence, it might also change the 
history of angular momentum transfer process which plays a 
major role in hierarchical scenarios (e.g. Navarro et al. 2004; 
D08). 

We study a set of intermediate-resolution cosmological 
simulations where different baryonic structures have been 
formed from the same initial condition (IC), through the 
modification of the physics of baryons (Scannapieco et al. 
2008, S08). These experiments are taken as toy models to 
analyse how the dark matter evolves when baryons are as- 
sembled in a different fashion. In this letter, we focus our 
discussion on four experiments which better illustrate the 
main results (for further details see Pedrosa et al. in prepa- 
ration) . 



2 SIMULATIONS AND ANALYSIS 

We analysed four simulations of a « 10 12 Mq type galaxy 
run with an extended version of the code GADGET-2 (Scan- 
napieco et al. 2005, 2006). Although all cosmological nu- 
merical simulations to date suffer from a number of lim- 
itations such as coarse-resolved star formation and inter- 
stellar medium physics, this extended GADGET-2 code has 
been designed to improve the respresentation of the ISM 
and SN feedback (S08). Our simulations are fully cosmolog- 
ical and include a multiphase model for the gas component, 
metal-dependent cooling, chemical enrichment and energy 
feedback by type II and la SN. The IC corresponds to a ^ 
10 12 M© halo extracted from a cosmological simulation and 
re-simulated with higher resolution. This halo was required 
to have no major mergers since z — 1. The simulations are 
consistent with a ACDM Universe with Q\ = 0.7, Q m = 0.3, 
Q b = 0.04, cr 8 = 0.9 and H = 100/i km s" 1 Mpc" 1 , 
with h = 0.7. The particle mass is 1.6 x 10 7 /i _1 M© and 
2.4 x 10 6 /i -1 M© for the dark matter and baryonic particles, 
respectively. The maximum gravitational softening used was 
e g = 0.8/i -1 kpc. At z = 0, the haloes are relaxed as indi- 
cated by a relaxation parameter of ~ 0.002 (Neto et al. 
2007). All analysed DM haloes have more than 120000 par- 
ticles within their virial radius (the total number of particles 
goes from w 170000 to « 320000 depending on the star for- 
mation histories). We analysed three simulations (NF, E-0.7 
and E-3) of the same IC but using different input param- 
eters for the SN feedback model. As a result, DM haloes 
finally host baryonic structures with different morphologies. 
In particular, simulation NF has been run without including 
SN feedback, while E-0.7 and E-3 include feedback assum- 
ing an energy per SN of 0.7 x 10 51 and 3 x 10 51 erg, re- 
spectively. For comparison, we also performed a pure grav- 



Table 1. Main characteristics of DM haloes and the main cen- 
tral galaxy. We show the total stellar mass M s of the central 
galaxy, the total-mass-to-stellar mass ratio Mt/M s , the n and 
r_2 Einasto parameters and the inner logarithm slope 7i nn er- All 
masses are evaluated at two times the radius that enclosed 83% 
of the baryons in the central regions. Bootstrap errors for n and 
r_2 are shown within parenthesis. 



Run 


10 10 h _1 M 


M t /M s 


/ ; 


r-2 
kpc h~ 1 


7inner 


NF 


15.9 


4.1 


5.973 (1) 


16.90 (1) 


1.24 


E-0.7 


7.5 


5.1 


6.887 (4) 


15.36 (1) 


1.28 


E-3 


1.3 


37.2 


5.585 (1) 


21.89 (1) 


1.15 


DM-only 






5.239 (3) 


24.06 (1) 


1.08 


DNF 


0.018 


2.4 


7.164 (5) 


1.10 (1) 


1.30 


DE-0.7 


le-4 


1486 


5.445 (4) 


2.17 (1) 


1.00 


D-only 






5.222 (3) 


2.41 (1) 


1.00 



itational run (DM-only) of the same IC (particle mass of 
1.84 x Mq). 

A full discussion of the effects of SN feedback on the 
baryonic dynamics and star formation histories can be found 
in S08. Here, we only summarize the main characteristics 
of the galaxies which are relevant to our discussion. At 
z — 0, the galaxy formed in NF is dominated by a extended 
spheroid, with most of its stars formed at z > 2. E-0.7 has 
been able to produce a galaxy with an important disc com- 
ponent as a result of the regulation of star formation by SN 
feedback. This system has a half mass radius of 5.72 h~ x kpc 
and a mass disc-to-spheroid ratio of 0.82. In E-3, the large 
amount of energy assumed per SN triggers violent outflows 
which expel a significant amount of the gas content of the 
main galaxy, producing a small thick stellar disc. Table [1] 
lists the final stellar masses of the galaxies formed in these 
simulations. 

We constructed the DM profiles for the simulated galax- 
ies at z — 0, after cleaning them from substructures. We fit- 
ted the NFW, Jaffe (1983) and Einasto (1965) expressions to 
the DM profiles, between three times the gravitational soft- 
ening and the virial radius, finding that the Einasto model 
provides the best fit in all cases. The Einasto formula has 
three free parameters, n, r_2 and p_2, which indicate the 
sharpness of the profiles, and the radius and density where 
their logarithmic slope takes the isothermal value. The fit- 
ting values obtained for the different simulations are shown 
in Table [1] We estimated bootstrap errors for n and r_2 by 
fitting the Einasto's formula to 100 randomly-generated re- 
alizations of the DM profiles and by estimating the standard 
dispersion over the generated set of parameters. 

When baryons are considered, the shape of the DM 
profiles in the central regions changes significantly in com- 
parison to the DM-only case, as it is shown in Fig. [T] The 
profile of the DM-only run is sharper (i.e. smaller n value) 
than those of haloes with baryons, indicating the increase 
in the DM concentration in the latter cases. In E-3, where 
an extreme (and perhaps unrealistic) value of 3 x 10 51 erg 
was assumed, we obtained the less concentrated DM pro- 
file among the cases with baryons. In fact its DM profile 
is weakly more concentrated than the DM-only one. As ex- 
pected, the galaxy in E-3 is also the most DM dominated 
one in the central region, as indicated by the total to stellar 
mass ratio Mt/M 8 (Table [TJ). In this simulation, most of the 
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Figure 1. Dark matter profiles for the NF (solid line), E-0.7 
(short-dashed line), E-3 (dotted-dashed line) and DM-only (thick 
long-dashed line) experiments. The inner bin corresponds to three 
gravitational softenings. The errorbars have been estimated by 
boostrap resampling technique. 



gas has been blown away and only a small fraction of stars 
has been formed in a bursty fashion. 

The comparison between the E-0.7 and the NF profiles 
shows that the DM distribution in E-0.7 is slightly more 
concentrated than in NF, although it hosts a galaxy a fac- 
tor of two less massive than the later. This finding suggests 
that the total amount of baryons collected within the cen- 
tral region of a halo is not the only relevant factor affecting 
the response of the DM to the presence of baryons (this be- 
haviour is common to the full S08 simulation sample as it 
will be discussed in Pedrosa et al. in preparation). The NF 
and E-0.7 runs have produced very different galaxies as we 
have already explained. In Fig|2] we show their age-radial 
position maps where the star formation histories and spa- 
tial distribution of the stars can be clearly appreciated. The 
galaxy in the NF run is dominated by old stars, determin- 
ing an extended spheroid with 78% of the final stellar mass 
older than 10 Gyr. The system in E-0.7 has a compact old 
spheroid and an important disc component populated by 
younger stars. This disc survives the interaction with satel- 
lites at lower redshifts (S08). The different star formation 
histories and morphologies in NF and E-07 are the result of 
the action of the feedback which, in the E-0.7 regulates the 
transformation of gas into stars preventing the formation of 
an early, extended spheroid and assuring the existence of gas 
to form a disc later on. The SN feedback has also affected 
the formation of stars in the satellite systems accreted by the 
main system, so that in the NF run they are more massive. 

We compare the residuals of the circular velocities esti- 
mated by applying the AC method of B86 to the DM-only 
run taking into account the baryonic distributions of the 
corresponding hydrodynamical runs (thin lines) and those 
obtained directly from the simulations (thick lines) as shown 
in Fig. [3] All residuals are estimated with respect to the cir- 
cular velocity of the DM-only run. We found that the B86 
model not only overpredicts the level of concentration but 
also changes the shape of the DM distribution when com- 
pared to the DM profiles obtained from the cosmological 
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Figure 2. Age-radial distance maps of the stars in the NF (upper) 
and E-0.7 (lower) experiments. This figure shows the inside-out 
formation of the disc in E-0.7 and the outside-in formation of the 
spheroid in NF. 



runs. While in the B86 method, the effects are principally 
determined by the amount of baryons accumulated at the 
centre, the simulated profiles are the result of the joint evo- 
lution of baryons and DM. More sophisticated AC meth- 
ods could account for these differences as suggested by the 
results of Gnedin et al. (2004) and Sellwood & McGaugh 
(2005). We will explore this possibility in a forthcoming pa- 
per. 

As can be seen in Fig. [3] while the simulations yield the 
largest residuals for the E-0.7 run which has the important 
disc component, the B86 model predicts the strongest effect 
for the NF run because it has the largest baryonic mass sys- 
tem, almost twice more massive than in E07. In the case 
of E-3 similar level of residuals are obtained. In the very 
central regions (within one e g and three e g ) our simula- 
tions predict an expansion of the DM profiles when baryons 
are present which can be related to the evolution of the sys- 
tems as discussed below. We also found this behaviour in the 
higher numerical resolution haloes of the Aquarius projects 
(Springel et al. 2008b). 

These results suggest a connection between the DM evo- 
lution and the history of formation of the baryonic struc- 
tures. To study this point, we estimated the concentration 
parameter proposed by Alam, Bullock & Weinberg (2002), 
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Figure 3. Residuals from the DM circular velocity measured in 
the simulations (thick lines) and from the AC prescription (thin 
lines) applied to the pure dynamical run considering the baryonic 
distribution in the NF (solid lines), E-07 (dashed lines) and the 
E-3 (dotted-dashed lines) simulations. The arrows denotes one 
(short) and three (long) gravitational softenings 
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Figure 4. Central halo mass concentration A v / 2 as function 
of redshift for the NF (solid line), E-0.7 (short-dashed line), E- 
0.3 (dotted-dashed line) and DM-only (long-dashed line) experi- 
ments. The solid line is the expected growth due to the expansion 
of the Universe. 

This parameter has the advantage of being indepen- 
dent of a specific density profile and, as it is an integrated 
quantity, it can be estimated more robustly at any time. In 
Fig. |4]we show A v / 2 as a function of redshift for our set of 

1 This parameter measures the mean DM density normalized to 
the cosmic closure density within the radius at which the circular 
rotational speed due to the DM alone rises to half its maximum 
value. 



simulations. All haloes increase their central concentration 
as they grow with time. The solid line shows the expected 
relation (i.e. hereafter critical relation) due to the expan- 
sion of the Universe alone. The DM-only run has the lowest 
DM concentration, as expected, at all times, although it flat- 
tens between z « 1 and z ~ 1.8, coinciding with the close 
approaching of two satellites. The transfer of angular mo- 
mentum from the satellites to the dark matter can explain 
its subsequent expansion (Pedrosa et al. in preparation). In 
E-3, A v /2 grows following the slope of the critical relation 
even closer than in the DM-only run since the entrances of 
satellites cause weaker effects on the DM distribution be- 
cause they are less massive due to the strong action of SN 
feedback. An interesting behaviour is determined by the NF 
and E-0.7 runs. Both haloes have similar concentrations un- 
til z ~ 1.5, and from there on, the E-0.7 halo has systemat- 
ically higher values of A v / 2 than the NF one. The analysis 
of the baryonic dynamics of these tests shows that by this 
redshift the disc component in E-0.7 starts forming. The in- 
coming satellites are less massive and the gas is able to settle 
onto a disc smoothly. This and the fact that the stellar disc 
is stable produce a smaller effective transfer of angular mo- 
mentum to the dark matter which, at the end, is slightly 
more concentrated than in the NF. 

The shape in the central region of the DM profiles can 
be also quantified by calculating the logarithm slope j(r) = 
— dlnp/dln r at the innermost resolved radius. We found 
7inner ~ 1 for the DM-only run (Table [T]), in agreement with 
recent estimations of Navarro et al. (2008) for galaxy-sized 
DM haloes in the Aquarius Project. The presence of baryons 
increases 7mner to « 1.25 as it can be seen in Table [1] 

The shape of the dark matter profiles in dwarf galaxies 
is one of the main caveats of the CDM cosmology. And it is 
also a difficult problem to tackle by numerical codes due to 
limited numerical resolution and oversimplified modelling of 
baryonic physics. To shed light on this point, we studied a 
scaled-down version of our « 10 12 M© type galaxy designed 
to get a system of 1O 9 M0. The maximum physical softening 
used in this run is 80 pc h -1 . This experiment (DE-0.7) has 
been run with the same cosmic and feedback parameters as 
E-0.7. However, the effects of the SN feedback on baryons 
are much stronger because of its shallower potential well as 
explained in detail by S08. We compared the DM profiles 
of this DE-0.7 run with those obtained in a no feedback 
(DNF) case and in a pure dynamical run (D-only) of the 
same scaled- down galaxy. 

For this small galaxy, SN feedback strongly prevents 
the star formation activity compared with that of the DNF 
case. In this shallow potential well, violent winds are able to 
blow out an important fraction of the baryons, inducing the 
system to form stars in a series of very weak starbursts after 
an early important one. As a consequence the DM profile 
is more similar to that of the pure dynamical run (Fig. [5}. 
The inner logarithm slope of the DM profile in the DE-07 
and DNF are 7mner = 1.0 and 7mner = 1.3, respectively. Con- 
versely in the DNF case the central accumulation of baryons 
induced a higher concentration of DM compared to the D- 
only run. As shown in Table 1, the shape parameter of the 
Einasto model goes from n — 5.2 in the D-only to n — 5.4 
and n = 7.2 in the DE-07 and DNF, respectively. We found 
that the evolution of A v / 2 with time of the DM profile in 
the DE-0.7 is very similar to that of the pure dynamical 
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Figure 5. Dark matter density profiles for the scaled-down sys- 
tem of 10 9 /i _1 M virial mass, for the DNF (solid lines), DE-0.7 
(short-dashed lines) and the D-only (long-dashed lines) experi- 
ments. The inner bin corresponds to three gravitational soften- 
ings. The inset plot shows the corresponding A v /% concentration 
parameters as a function of redshift. 



run in agreement with previous results (Read, Pontzen and 
Viel 2006). This could be caused by the too strong SN feed- 
back effect which blows away most of the baryons or by 
deficient numerical resolution. Indeed, by using very high 
numerical resolution simulations and considering molecular 
cooling, Mashchenko, Wadsley and Couchmanet (2008) de- 
tected strong differences in the DM profiles at z = 5. 



3 CONCLUSIONS 

We found that, in all our 10 12 haloes, including the baryons 
resulted in a more concentrated DM density distribution 
compared to the equivalent pure DM simulation. However, 
using the total baryonic mass that collects at the centre as 
input to the AC model of B86, we were not able to reproduce 
our final DM density distributions. Our results indicate that 
the way baryons are assembled into a galaxy affect the DM 
evolution so that its final distribution also depends on the 
morphology of the baryonic system. An efficient star forma- 
tion activity produces satellites with larger stellar masses 
which are able to survive further in the halo. This could 
cause a larger angular momentum transfer (D08) which, on 
its turn, might explain the expansion of the DM central 
concentration. The treatment of SN feedback prevents this 
early overproduction of stars. Hence the gas is allowed to get 
to the centre before being transformed into stars forming a 
more compact stellar spheroid at early times. At later times, 
the remanent gas and the disrupted satellites are able to set- 
tle onto more important disc structures. When SN feedback 
is very efficient and blows away significant masses of gas, 
then the DM profile is less affected because of both the low 
central baryonic concentration and the easier dissolution of 
smaller satellites that get into the halo. We found that this 
behaviour can be linked to a larger DM concentration pre- 
sumedly due to a weaker angular momentum transfer from 
the satellites to the dark matter. Our findings confirm pre- 
vious lower resolution results and suggest new clues for the 



understanding of galaxy formation, which will be explored 
in a forthcoming paper. 
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